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Propagation of an extraordinary electromagnetic wave in the vicinity of electron 
cyclotron resonance surface in an open linear trap is studied analytically, taking into 
account inhomogeneity of the magnetic field in paraxial approximation. Ray trajec- 
tories are derived from a reduced dispersion equation that makes it possible to avoid 
the difficulty associated with a transition from large propagation angles to the case 
of strictly longitudinal propagation. Our approach is based on the theory, originally 
developed by the Zvonkov and Timofeev [1], who used the paraxial approximation 
for the magnetic field strength, but did not consider the slope of the magnetic field 
lines, which led to considerable error, as has been recently noted by Gospodchikov 
and Smolyakova 2j. We have found ray trajectories in analytic form and demon- 
strated that the inhomogeneity of both the magnetic field strength and the field 
direction can qualitatively change the picture of wave propagation and significantly 
affect the efficiency of electron cyclotron heating of a plasma in a linear magnetic 
trap. Analysis of the ray trajectories has revealed a criterion for the resonance point 
on the axis of the trap to be an attractor for the ray trajectories. It is also shown 
that a family of ray trajectories can still reach the resonance point on the axis if the 
latter generally repels the ray trajectories. 

As an example, results of general theory are applied to the electron cyclotron 
resonance heating experiment which is under^ preparation on the Gas Dynamic Trap 
in the Budker Institute of Nuclear Physics 



PACS numbers: 52.50. Sw, 52.25.Xz, 52.40.Db 



2 



I. INTRODUCTION 

The behavior of electromagnetic waves in plasmas in the range of the electron cyclotron 
frequency remains the subject of research for a long time. The main results of these studies 



are well known — see, for example. 



7|, and the interest in the problem is supported at high 
level due to important applications such as high-frequency heating and plasma diagnostics 
in magnetic traps. 

The theory of propagation of electromagnetic waves in a plasma with a one-dimensional 
inhomogeneity is developed in detail. In particular, in a plane-layered system, spatial res- 
onance is a singular point of the wave equation [7]. In terms of ray trajectories, it means 
that the movement of a wave packet slows down along the inhomogeneity as the packet 
approaches the point (plane) of spatial resonance. This leads to the accumulation of wave 
energy in the vicinity of the resonance point and the appearance of singularities in the 
distribution of electromagnetic fields. 

In the case of two-dimensional inhomogeneity, determining the position of the resonance 
surface is more complicated. According to the pioneering paper it can be done by 
analyzing the trajectories of wave packets (ray traces). 

Due to the short wavelength in the range of frequencies of electron cyclotron resonance, 
the wave propagation in laboratory magnetic traps is described with reasonable accuracy 
by the approximation of geometric optics, and the investigation of the ray trajectories can 
provide both qualitative and quantitative assessment of the possibility of plasma heating in 
the trap. 

In this paper, propagation of an extraordinary electromagnetic wave in the vicinity of 
electron cyclotron resonance surface in an open linear trap is studied analytically with tak- 
ing into account the heterogeneity of the magnetic field in paraxial approximation. Wave 
trajectories are described by solving a reduced dispersion relation that makes it possible 
to avoid difficulties associated with a transition from large propagation angles to the case 
of strictly longitudinal propagation. Our approach is based on the theory, initially devel- 
oped by Zvonkov and Timofeev jl], who used paraxial approximation for the magnetic field 
strength. Unfortunately, they ignored the inclination the magnetic field lines, which resulted 
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in significant error as was recently noted by Gospodchikov and Smolyakova 2\. We refine 
the theory of Ref. |l] in order to take into account the effect specified in Ref. 2|. The refined 
theory makes it possible to find ray trajectories in analytic form. 

As an example, we apply the results of general theory to the electron cyclotron resonance 
heating (ECRH) experiment at the Gas Dynamic Trap (GDT) which is under preparation 
in the Budker Institute in Novosibirsk [3]. The key physical issue of the GDT magnetic field 
topology is that conventional ECRH geometries are not accessible. The solution, proposed 



and numerically investigated in [3], is based on a peculiar effect of radiation trapping in 
inhomogeneous magnetized plasma. Under specific conditions, oblique launch of gyrotron 
radiation results in generation of right-hand-polarized (R) electromagnetic waves propagat- 
ing with high longitudinal refractive index N\\ in the vicinity of the cyclotron resonance 
layer, which leads to effective single-pass absorption of the injected microwave power. The 
theory developed in the present paper is related to the final stage ECRH when the trapped 
electromagnetic wave approaches the resonant surface, where the wave frequency u is equal 
to the local electron cyclotron frequency u e (r). 

This paper is organized in the following way. To begin with, in Sec. [IT] we describe 
the magnetic field in a linear trap in the paraxial approximation. In Sec. II I II we derive 
the dispersion equation for electromagnetic waves near the electron cyclotron frequency. 
In Sec. IIV[ we find parametric formulas for the ray trajectories and derive a criterion for 
the resonance point on the plasma axis to attract the ray trajectories. In Sec. |V] the ray 
trajectories are analysed and systemized. Finally, in Sec. IVIl the effect of finite plasma 
pressure on the electron cyclotron resonance heating is discussed. A short summary is given 
in Sec. Em 

II. MAGNETIC FIELD IN PARAXIAL APPROXIMATION 

Neglecting the the finite (3 effects (where /3 is the ratio of plasma pressure to the pressure 
of magnetic field), we assume that the magnetic field in the plasma is potential, that is 



B = W, 



(1) 



4 



where the potential ip obeys the Laplace equation V 2 ip = 0. For an axisymmetric system in 
cylindrical system of coordinates the latter has the form 



Putting here 



yields a recursive formula 



r dr \ dr ) ~^ dz 2 ^ ^ 



i> = Y,Mz)r n (3) 



< + in + 2fA n+2 = 

for the functions A n (z). It allows to express all the functions A n (z) either through A (z) 
or Ai(z). Since the radial component B r (r, z) of B is zero at r = due to axial symmetry, 
the coefficient Ai(z) should be zero, and the coefficient A (z) can be expressed through the 
magnetic field Bq(z) = B z (0,z) on the system axis. Paraxial approximation means that 
high-ra terms in ([3]) are small and can be dropped. Keeping first two non-zero terms of the 
expansion, we obtain 



V>(r) 



J dzB (z) - jB' (z), (4) 



where the prime stands for the derivative over z. Consequently, 

B r (r) = ^ = - r -B' (z), B z (r) = ^ = B (z) - jB'^z). (5) 



Following Ref. [lj], we expand the modulus of the magnetic field strength 

B(t) = v^WTBl « B (z) - ^ {b'i - \$L\ (6) 



4 V 2 B Q f 

about the point z = z s of the electron cyclotron resonance on the system axis to put it in 
the form 

B(v) = B (z s ) (l + - , (7) 

where 

1 B' (z a ) 1 _ l Bg(z B ) l B>£{z s ) 

L,| B (z s y L\ 4B (z s ) 8% s )' U 
As an example of general theory implementation, we consider ECRH in GDT at the 
frequency u/2n = 54.5 GHz. We consider two magnetic field configurations, which differ by 
the maximum of the field strength B max at the magnetic mirrors. Low field configuration 
is characterized by -B max = 117 kGs, and strong field configuration has B m3bX = 145 kGs. 
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FIG. 1: On-axis distribution of the magnetic field strength Bq(z) in the GDT device for two con- 
figurations with maximum magnetic field 117 kGs and 145 kGs; the dashed line shows the magnetic 
field strength corresponding to the electron cyclotron resonance at the frequency 54.5 GHz. 



oj/2-k (GHz) 


S max (kGs) 


z s (cm) 


Ln (cm) 


L\ (cm 2 ) 




54.5 


116.9 


327.5 


29.8 


770.7 


1.16 


54.5 


145.4 


306.6 


136.3 


-3768.0 


-4.93 



TABLE I: The characteristics of the electron cyclotron resonances in GDT. 

Fig. [T] shows plots of the magnetic field profile B (z) in the west half of GDT for these two 
configurations, and the horizontal line indicate the resonant magnitude of the magnetic field 
B res = 19.5 kGs. Resonance points z s are located at the intersection of this line with the 
plots of B (z). Computed values of the parameters z s , L\\, L\ are given in Table [H 

III. RAY HAMILTONIAN 

To derive the ray hamiltonian, we shall use the helical (spiral) coordinates 

± dx±idy I, 

as = -= — , as" = az. 

v2 

The metric tensor can be found from the identity 

ds 2 = 2 ds + ds~ + ds" ds" = g ik ds l ds k 
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and has the form 

f Q 1 o\ 

9ik= 10 

V 001 / 

which immediately yields the rule ds± = ds T , dsy = ds" for the covariant components 
dsj = gikds k . Components of other vectors are defined similarly. In particular, helical 
components of the electric field E and normalized wave vector N = ck/u (refractive index) 
are 



E ± = E X ± iE y 



V2 



E. 



E" = E, 



N 



± 



N x ±iN v 

y _ N 



V2 

Contravariant components of the permittivity tensor e a/3 are computed through those in the 

cartesian coordinates by the standard rules. For example 

++ ds + ds + 1 . . 

8 = ~dx^~dx^ £i ^ = 2 x ~ £yy ^ l£xy ~^ l£yx ' ' 
In the cold plasma limit, the diagonal elements e ++ and e are zero whereas the off-diagonal 
components are not: 



+ _ ds + ds 1 



_ + ds ds + 1 . . 

£ = ~d^~dx^ Sii = 2 x ~^ £yy ~^ %£xy ~ %£yx ' = 
All other components except for = e zz = £y are also zero so that 

^0 e + 0^ 

en 



T af3 



where 



UJ 



ps 



UJ. 



pe 



^ UJ 2 ~ UJ 2 



ui{ui =f Q s ) w(u;±u; e )' 
summation goes over the plasma species s, oo ps = \jAne 2 n s /m s , Vl s = e s B/m s c, and co e - 
— Q e = \e\B/m e c. Mixed permittivity tensor is diagonal in the cold plasma. In particular, 

U + 0^ 



(9) 



e_ 
e 
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In the approximation of geometric options, the wave equation reads 

{e% - N 2 g a p + N a Np) E p = 0. 

It leads to the dispersion equation {?J 

N + N + N~N- N\\N\\ , , 

and the polarization of the eigenmodes is described by the vector 

r /v+ n- n\\ \ 

E^*-,*« } oc{_^,_^,_ - 2 ). (11) 

Note that N + N + = N~N- = |(iV| + A" 2 ) = |iVj_, and E + and E~ are the amplitudes of 
the left and right circular wave components that rotate in the direction of the electron and 
ion cyclotron gyration, respectively. 

Below we restrict ourselves to the case to ~ ui e , where e_ — > oo. 

Assuming that A^ 2 <C we can drop the third summand in Eq. (flOl) . Then, for trans- 
verse propagation, N\\ — > 0, we find two eigenmodes. The ordinary mode is polarized along 
the magnetic field: 

N 2 ~e {l , E~ {0,0,1}. 

It is obtained from the dispersion equation (1101) by evaluating the 0/0 uncertainty in the 
last term. The extraordinary mode is circularly left polarized (and gyrates in the ion side): 

N 2 ~2e + , E~ {1,0,0}. 

In the case of longitudinal propagation, 9 = N±/N — > 0, the ordinary node can still be 
found using the assumption iV 2 <C It has the left (ion) circular polarization: 



iV 2 ~ e. 



2en 



E~ {1,0,0}. 



Dispersion and polarization of the extraordinary wave undergo significant changes as the 
angle 6 approaches zero. First of all it is worth noting that iV 2 = e_ is the exact solution at 
9 = and, hence, iV 2 — > oo as £_ — > oo so that the assumption iV 2 <C breaks. Proceeding 
with greater accuracy as demonstrated in Ref. [1|, we put 

iV+AL N 2 /2 NWNh Nl-ea 

e+ — N 2 N 2 ' en — A^ 2 N 2 
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into Eq. (ITUj) to obtain the reduced dispersion equation 



q|iV 2 + \e-Nl -£||£_ = 0. 



(12) 



Its solution 



N 2 



£||£_ 




shows that the dispersion law abruptly changes at 




since N 2 



~ £_ 



for < 0, whereas N 2 ~ 2^ N /6> 2 and ~ 2er M for 9 X < < tt/2. 

In the range of frequencies u; > w e , where < 0, a plasma resonance, characterized by 



and the plasma is transparent only for 9 > 9 CI . It means that a packet of extraordinary wave 
can hardly access the surface of electron cyclotron resonance, where u = u e , from the side 
of low magnetic field. This fact will be confirmed by direct calculations in Sec. IVl 

For 9 ^> 9[, the wave field has almost complete longitudinal polarization (En E + ^> 
E-) since 



and the left circular component is small (E + <C E~,E^). The polarization drastically 
rearranges at the angle 9 ~ 9f. In the interval 9f <C 9 -C 9i, the field still has almost 
complete longitudinal polarization ;§> E- ^> E + ) as in the case 9 ^> 9\. On the 
contrary, for 9 <C 9f, the wave field has almost complete right circular polarization that 
rotates in the direction of the electron gyration (E_ ^> Eh ^> E + ). 

It is important to note that Nj_ remains finite in the limit £_ — > oo. This allows us 
to substitute the first term in Eq. (|T2|) with £\\N 2 since £||iV 2 is small as compared to the 
second term. Second simplification is due the observation that near the electron cyclotron 
resonance, at u ~ u e , it is sufficient to keep only the last term in the expression (j^J) for e_. 
Then, Eq. (|12p takes the form 



iV — > oo, occurs at the critical angle 9 CT = ^—e\\/2e_ (we assume that e\\ > 0, i.e. u > u; pe ), 




In the range 9 <C 6i, the polarization is given by the vector 




D = euN; 



2 u(u e (r) -u) 1 
11 < 2 




£|| = 0, 



(13) 
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a r 




FIG. 2: The effect of magnetic field lines inclination. 



where the factor 



[UJAr 



to 



z - z. 



u 



r 



(14) 



is computed using Eq. (J7J). 

In a weakly inhomogeneous magnetic field, Eq. f|T3|) defines slowly varying function D of 
the radius- vector r. Given the validity of geometric optics, the function D = D(k, r; u) can 
be considered as the ray Hamiltonian [5] that governs motion of a wave packets similar to 
the ordinary Hamiltonian that governs motion of a particle in mechanics. Appropriate set 
of equation 



dr 

dt 



3D 



dk 

dt 



3D 
dr 



(15) 



relates the radius-vector r and the wave-vector k of the packet with some "time" t which, 
however, has neither implicit meaning nor dimension of physical time. 

The transversal and longitudinal components, iVj_ and N\\ , of the vector N can be matched 
to the radial and axial projections, iV r and N z , as shown in Fig. [21 where ip ~ B r /B = 
—r/2L\\ is the inclination angle of the magnetic field lines: 



N± = cos ipN r — sin ipN z 



Nn = cos ipN z + sin^A^ r 



Since i/)<l and N± Nn we have approximately 



N ± w iV r - ipN z Nn w N z . 



(16) 
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To proceed further, we introduce a dimensionless parameter 

,2 



a — 



4£||u; 2 

and normalize coordinates on the parameter of length 



2aU 



'pe 



so that 

p = r/£, x = (z- z s )/£. 

Instead of the dimensionless vector N, we introduce the dimensionless vector of the wave 
momentum p = (u£/c)~N with the components 

oo£N r u£N z 



Pp = > Px = • 

c c 

With all these amendments being combined, the Hamiltonian (TX5]) takes the form 



where 



D = (x~ W 2 )p\ + ip P + otpp x ) 2 - v = 0, (17) 



T \ 2 , A T2 



V = 2a ' / 2 2 

c / u z c z 2e\\ 

In Ref. |l|] the term term app x in Eq. fTT7|) was lost since iV^ was set equal to N r . In what 
follows we assume that a > 0, i.e. w > cj pe . ECR heating is not effective in the case uj < u pe 
as will be discussed later. 



IV. STABILITY OF THE RESONANCE ZONE 

Following to Ref. [lj], we consider only those ray trajectories that lay in a plane passing 
through the axis z of the system. Appropriate ray trajectories for the Hamiltonian ffTTj) are 
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governed by the equations 

2( Pp + app x ) 1 (18) 

2p J pp x -2ap x (p p + app x ), (19) 

2(x~ W 2 )p x + 2ap(Pp + apPx)> (20) 

-PL, (21) 



dp 


dD 


dr ~~ 


dp p 




dD 


dr 


dp 


dx 


dD 


dr 


d Px 


dp x 


dD 


dr 


dx 



where r is the new dimensionless time. The last equation is separated from others and has 
the solution p x = l/(r — C\), which depend on the arbitrary constant C\. The latter can 
be set to zero by appropriate choice of the origin of time; hence 

Px = 1/r- (22) 

Putting (122!) into Eqs. (ITS"]) and (jl"9l) . we note that their solution on the interval r > can 
be sought in the form 

p = A+t^ + A_r 7 ", (23) 
Pp = B+t^- 1 + B^t'— 1 . (24) 



It is then easy to find that 



7± = ~ (l ± x/1 - 8a + 16/x) , (25) 
5± = i(7±-2 Qt )A t . (26) 
The function x can be found from Eq. (|T7|) : 

x = mp 2 + ^ - ( t Pp + «p) 2 = 

vt 2 + (4/i - a) A_ A + t - ~ (2 7+ - 4a) A^r 2 ^+ - - (2 7 _ - 4a) AV 7 - . (27) 

The obtained solution is valid for /? < 1, % < 1, since it is based on the paraxial ap- 
proximation Second condition of validity reads p P /p x — Pp T <C 1; it guaranties that 
= Pp r + < 1. In terms of the time r the validity region is limited by the inequalities 
\A + T~t+ \ < 1 and |A_r 7 -| < 1. 

The functions p(r) and x(r) determine a ray trajectory in the parametric form. They 
depend on the two amplitudes A + and A_, which can be related to the initial coordinates of 
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FIG. 3: (Color online) Profile of the ratio L^/L 2 ^ in the GDT device for the low (solid line) 
and strong (thin line) magnetic field configurations. Attr actor region is located below the dashed 
line. Vertical solid and dashed arrow show position of the fundamental cyclotron resonance at the 
frequency 54.5 HGz for the low and strong field configurations, respectively. The low-field-resonance 
falls into the repeller region, and the strong-field-resonance is located within the attractor region. 

the ray po and Xo & t a given initial value of the parameter r . We assume that a particular 
trajectory starts at some r > and ends at r = 0. Note that ECRH occurs at r — > when 
p x — > oo (i.e. N\\ — > oo) and the extraordinary waves transforms to small-scale oscillations 
that are effectively absorbed by plasma in single pass of the injected wave. 

Absorption of the wave energy on intermediate parts of the ray trajectory is weak even 
if the trajectory crosses the resonance surface, where u = co e ; the latter is described by the 
equation 

X = (28) 

in dimensionless notations. 

It is readily seen that both p and x tend to zero at r — > 0, if the real parts of both r y + and 
7_ are positive, Re7± > 0. This occurs if and only if 1 — 8a + 16/x < 1 and means that any 
ray trajectory reaches the resonance point p = x — on the axis of the system. In this case 
the point p = x — turns out to be an attractor for the ray trajectories. In the opposite 
case the resonance point repels most of the trajectories (see below). 

For a > the condition 1 — 8a + 16p < 1 is equivalent to the inequality 



T 2 i 



(29) 



The criterion f l2"9~j) lays in-between those obtained in Ref. [lj], where it was equivalent to 
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L 2 L\ < 0, and in Ref. {2], where the inequality L 2 L\ < ~ was derived from approximate 
reasonings. 

As an example, Fig. |3] plots the ratio L 2 /L\ for the two magnetic field configurations in 
GDT referred to in table [H Attractor region (1291) is located below the dashed line in Fig. [3l 
and the resonances are indicated by vertical arrows. It is remarkable that the two plots (thick 
and thin curves for the low and strong field configurations, respectively) almost coincide 
whereas the positions of the resonant point z s differ significantly for the two configurations. 
Solid arrow stands for the low magnetic field configuration and indicates the resonance 
point z s = 327.5 cm located in the repeller region. Dashed arrow indicates the position of 
the resonant point z s = 306.6 cm for the strong magnetic field configuration; it is located in 
the attractor region. 

Despite the fact that the resonance point on the axis for the low field configuration does 
not satisfy the condition (I211 . it can be reached by appropriately tuned ray trajectories as 
is shown in the next section. 



V. RAY TRAJECTORIES 



Depending on the parameter ( = 1 — 8a + 16p, which is included in the definition (125]) of 
the exponents r y±, ray trajectories can be divided into three types, considered sequentially 
in the following three subsections. 



A. The case < ( < 1 

If < ( < 1, the exponents r y± are real and obey the inequalities < 7_ < 7+ < 1. It 
means that the on-axis resonant point p = x — serves as attractor for the ray trajectories. 

To begin with, we consider a special case A + = 0. Then, it is possible to exclude the 
parameter r and express the function \ through p. The result has the form 

X = v P 2h --\(l--2a)p 2 , (30) 

where v = u/A 2 /" 1 ' . The parameter v is usually big and positive so that a typical trajectory 
approaches the resonant point at the axis from the side of strong magnetic field as shown in 
Fig. |H This observation is in agreement with the fact, noted in Sec. IHH that low magnetic 
field side is not transparent for the extraordinary waves. 
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FIG. 4: (Color online) Ray trajectories for A+ = 0, < ( < 1, a = 0.3, p = 0.1, v = 30, 7+ = 0.72, 
7_ = 0.28. Dashed line indicates the resonance surface uj = uj e . Ray trajectories approach the 
resonant point p = \ = from the side of strong magnetic field. 

P 



-0.001 




FIG. 5: (Color online) Ray trajectories for A + = 0, < C < 1, a = 0.125, p = 0.01, v = 10, 
7 + = 0.7, 7_ = 0.3. Dashed line indicates the resonance surface uj = uj e . Ray trajectories approach 
the resonant point p = x = from the side of low magnetic field but starts on the strong field side. 

Nevertheless there is a range of parameters, where the trajectory passes on the low field 
side just before accessing the resonant point as shown in Fig. Indeed, since 7_ < 1 the 
second term in Eq. (T5TTT) dominates in the limit p — > 0, where the trajectory is described by 
the parabola 

X = ~ (l - 4a - v 7 ! - 8a + 16//) p 2 . (31) 

Its branches are directed to the side of positive \ except for the narrow region of the param- 
eters a/2 — 1/16 < n < a 2 that exists if a < 7. In this narrow region the ray trajectories 
approach the resonant point p = \ — at the axis from the side of negative x (i- e - l° w 
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FIG. 6: (Color online) Ray trajectories for £ > 1, A- ^ 0, a = 0.5, p = 0.9, v = 5, 7+ = 1.17, 
7_ = —0.17. Dashed line indicates the resonance surface uj = u e . 

magnetic field) as shown in Fig. [5j 

If both coefficients A + and A_ are not zeros, qualitative picture remains the same and 
all the ray trajectories approach the resonant point p = \ = as they do in Figs. H] and 
It means that wave packages deposit their energy mainly in the plasma core. 

B. The case ( > 1 

If ( > 1, both 7+ and 7_ are real but now 7+ > 1 whereas 7_ < and the term A-T 1 - 
in Eq. fl23l) dominates at r — > except for the case where A_ = 0. In other words, all 
trajectories with non-zero coefficient A_ are repelled from the resonant point p = x — as 
shown in Fig. [6] so that the resonant point cannot be reached by most of the ray trajectories. 
Nevertheless, special ray trajectories with A_ = can still reach the on-axis resonant point 
and deposit there the wave energy as shown in Fig. [71 These trajectories are described by 
the equation 

X = vp 2h+ - -/ (7+ - 2a) , (32) 

where v = v/A 2 l~ 1+ . In the limit p — > 0, the first term in Eq. (|32i) dominates. Since it 
is positive, these special trajectories always approach the resonant point from the side of 
strong magnetic field. 

Another point to note for the case ( > 1 is where the repelled trajectories deposit the 
wave energy. As it was noted in Sec. \IV\ the energy is deposited into plasma at r — > when 
p x — > 00. In the case A_ 7^ 0, the radial coordinate p of a repelled trajectory formally tends 
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FIG. 7: (Color online) Ray trajectories for ( > 1, = 0, a = 0.1, = 0.8, v = 10, 7+ = 2.30, 
7_ = —1.30. Dashed line indicates the resonance surface uj = u: e . 

to infinity at r — > 0. It might occur that p x will reach sufficiently large values before the 
ray will leave the plasma and some part of the wave energy will be deposited at the plasma 
periphery. As can be seen from Fig. [6j in such a case the wave energy is deposited in the 
vicinity of the resonant surface f[28]) . 

The case of ( > 1 is realized in GDT for the low magnetic field configuration, where a 
varies from 0.02 to 0.5, \x from 0.015 to 1.1, and £ from 1.1 to 15, assuming that the plasma 
density in the heating zone varies from 0.5 x 10 13 to 2.5 x 10 13 cm -3 . 



C. The case ( < 

In the case ( < 0, the exponents 7+ and 7_ are complex conjugate to each other. The 
coefficients A + and should also compose a complex conjugate pair in order for p(r) to 
be real function. The latter can be rewritten as 

p(r) = A x t 112 cos (j v/l-8a + 16/x In A + A 2 r 1/2 sin (| y/l - 8a + 16// In A , (33) 

where A\ and A2 are the new amplitudes. The functions x( T ) an d Pp( T ) can then be found 
from Eqs. ffTTj) and ([TBI) , respectively. They describe the oscillating trajectory, which ends 
in the resonant point on the plasma axis. Consequently, this point turns out to be attractor 
for the ray trajectories as in the case < £ < 1, and the wave energy is mainly deposited 
near the attractor. 

Fig. [8J shows an example of such trajectories. They can approach the resonance point 
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FIG. 8: (Color online) Ray trajectories for ( < 0, a = 2, [i = 0.5, v = 10, Re7± = 0.5, Inry-i- = 
±1.32. Dashed line indicates the resonance surface uj = uj e . 

mainly from the side of strong magnetic field while some trajectories transfer on the weak 
magnetic field side just before entering the on-axis resonance point. 

In the GDT facility, the case £ < takes place for the strong field configuration as can be 
seen from table HI Numerical analysis of the ray trajectories for the strong field configuration 
performed in Ref. jsj] confirms the above conclusion that this configuration provides very 
effective ECRH. 



D. Accessibility of the on-axis resonance 

Suppose that Q < 1 so that an arbitrary ray trajectory ends near the on-axis resonance 
point p = x — 0- K does not mean however that the point can be accessed from anywhere. 
Indeed, if we take an arbitrary point (p , xo) and arbitrary time r , we conclude from Eq. (fT7j) 
that they are related by the equation 

Xo = PPo + UT o - ( t oPpo + «Po) 2 , (34) 

where 9q = Top p o + apo stands for initial value of the angle 9 w N±/N, which should be real 
for propagating waves as well as po and xo- Assuming that a > and, hence, v > 0, we 
conclude from Eq. (j34|) that for any 

Xo > PPo - #o ( 35 ) 

there exist such a time To and a ray trajectory starting from a given point (po, xo) at a given 
angle 8q, which reach the resonance point. Note that a trajectory with the angle — 6q of 
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opposite sign also reaches the resonant point from the same start point (po, Xo) f° r the same 
time t . 

The inequality (135]) implies that a trajectory can for sure access the resonant point if it 
starts on the strong field side x > PP 2 °f t ne resonant surface x = PP 2 ■ As it is explained 
in Sec. 1111} the low field side is not transparent for the extraordinary waves if a > (i.e. 
u > u pe ). 

In overdense plasma, where a < (u < uj pe ), the parameter v becomes negative, and 
instead of (ESJ) we formally obtain the opposite inequality Xo < PPo ~ $o> "which means that 
strong field side becomes opaque for the extraordinary waves, and the resonance point can 
be accessed mostly from the low field side. However the plasma periphery (where the plasma 
density is low and u > u pe ) will be opaque for the extraordinary waves in this case so that 
effective ECRH of overdense plasma is not possible as it is well known (see, e.g. |3|). 



VI. EFFECT OF FINITE /3 

So far we have ignored the effects that appear due to distortion of the magnetic field by 
the plasma pressure. Let now show that these effects are not important. 

Equilibrium of plasma cylinder assumes that the sum of plasma pressure and the pressure 
of the magnetic field is constant across the plasma, i.e. 

R2 R2 

It follows from Eq. (156]) that small plasma pressure changes the magnetic field inside the 
plasma by small amount 

SB = -f. (37) 
Assuming that the plasma pressure has parabolic radial profile, 

P = Po [l-- 2 

we conclude from Eq. ( |37|) that the parameter L\ should be substituted with 

Lf~ L\ 2a 2 ' {6 * } 

where (3 = 8t[p /Bq, and a denotes the plasma radius in the heating zone. As a result, the 
criterion ( 12"§|) takes the form 

Lm L 2 , / ar.2 
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It is seen from Eq. ( 13 9 p that the effect of plasma pressure is negligible if /3L^/2a 2 <C 1. In 
GDT, the factor (3L 2 ± /2a 2 does not exceed 0.01 and, hence, the effect of plasma pressure can 
indeed be neglected. 



VII. CONCLUSION 

In this paper we have refined the theory initially developed in Ref. jl| in order include 
the effect of magnetic field lines inclination as suggested in Ref. Q . We have found analytic 
form for the parametric presentation of the ray trajectories near the axis of axially-symmetric 
linear trap and derived the criterion ( I2"9~j) that guaranties for all the ray trajectories to reach 
the point of electron cyclotron resonant on the trap axis. Our criterion (129 j) lays in-between 
those obtained in Ref. jl|, where the effect of of magnetic field lines inclination was missed, 
and in Ref. (2j, where similar criterion was derived from approximate reasonings. 

We have found that topology of the ray trajectories qualitatively changes as the inhomo- 
geneity of both the magnetic field strength and the field direction varies, which significantly 
affects the efficiency of electron cyclotron heating of a plasma in a linear magnetic trap. 
In a case where the criterion f[2"9~j) is satisfied, all the ray trajectories end near the on-axis 
resonant point which turns out to be attractor. In the opposite case, the on-axis resonant 
point repels the most of the ray trajectories but, still, there is a family of the ray trajectories 
that ends at that point, which guaranties effective ECRH. 

As an example, we applied our theory to forthcoming ECRH experiment on the GDT 
facility JsJ. Our theory leads to the conclusion that the efficiency of ECRH is very sensitive 
to the position of the on-axis resonant point which should be located within the attractor 
zone in order for the plasma interior to be heated. We have shown that the strong magnetic 
field configuration which was numerically analyzed in Ref. |3(] meets this requirements. 
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